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AB8TRACT 

A  numrical  mthod  for  boundary  valua  problems  for  quasilinear  systama  of 
singularly  part ur bad  ordinary  differential  aquations  is  presented.  The  Method 
is  based  an  collocation  with  polynomial  splines.  The  stability  properties  of 
the  associated  difference  operator  are  examined  and  a  stepsize  algorithm  to 
achieve  a  certain  over-all  accuracy  is  developed.  The  number  of  gridpoints 
required  by  the  algorithm  is  estimated. 


AMS  (MOS)  Subject  Classifications!  34E99,  65L10 

Key  Wordst  Collocation  Methods,  Singular  Perturbations,  Boundary  Value 
Problems 

Work  Unit  Number  3  -  Numerical  Analysis 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 .  This 
material  is  based  upon  work  supported  by  the  Rational  Science  Foundation  under 
Grant  No.  MCS-7927062,  Mod.  2. 


SIGNIFICANCE  AND  EXPLANATION 


Many  high-order  discretisation  methods  for  ths  solution  of  two-point 
boundary  valus  problems  for  systems  of  ordinary  diffsrsntial  aquations  srs 
alrsady  circulated.  However,  these  asthods  can  behave  quite  poorly  in  case 
the  solution  has  large  derivatives,  unless  severe  restrictions  on  the  aesh  are 
imposed.  The  reason  for  these  restrictions  is  aainly  a  stability  problea. 

In  this  paper  a  strongly  A-stable  difference  aethod  based  on  polynoaial 
collocation  is  developed  for  a  class  of  quasilinear,  singularly  perturbed, 
two-point  boundary  value  problems.  Many  problems  of  practical  interest  are 
included  in  this  class,  for  instance  the  nonlinear  deforaation  of  thin  beams 
or  one-dimensional  models  of  carrier  transport  in  semiconductor  devices.  The 
aethod  combines  the  advantages  of  having  the  same  stability  properties  as  the 
lower  order  methods  which  are  used  already,  with  the  high  order  of  convergence 
of  collocation  methods.  It  is  shown  that  the  number  of  gridpoints  (and 
therefore  the  amount  of  computing  time  and  required  storage)  is  of  the  saae 
order  of  magnitude  as  the  one  required  for  solving  unperturbed  problems. 


The  responsibility  for  ths  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


ON  COLLOCATION  SCHXMXS  FOR  QUASILIHEAR 
SINGULARLY  FBRTURBCO  BOUNDARY  VALU1  FROBLBM8 


CHRISTIAN  RINGHOrtR 

1.  Introduction. 

Me  consider  quasilineer  systees  of  singularly  parturbad  O.D.B. '#  on  tha  intarral 
(0,1]  of  tha  fora 


(1.1) 


re<y.*> 

re<y.*> 

be(y,«) 


0  ,  be(y,a>  -  0 

ey1  (t,e)  -  f(t,y,s,c) 

*•  (t,e)  -  g(t,y,s,c) 
b(y(  0,  E),y(  1,0,  s(0,e),  s(1,C),C) 


y  and  a  ara  vectors  of  diaansion  n  and  a.  Tha  prime  danotaa  tha  darlvativa  with 
respact  to  t.  (1.1)  is  qoaailinaar  in  tha  following  sansai 


f(t,y,s,c)  “  f  (t,z)y  +  f  <t,y,z,e) 


(1.2) 


»>(y0.y, *  b1(«0,«1>(  °)  -*•  b2(y().y1,s0,*1,e) 


f.  and  b.  ara  aatricaa  of  appropriate  dimensions.  He  assume  that  tha  derivatives  r~ 

*y 

3bJ 

and  are  uniformly  of  order  e  for  y,s,yQ,y1  xQ,z1  in  any  bounded  domain  of  tha 

appropriate  real  spaces.  For  fj.bj  and  g  there  exist  asymptotic  expansions  in  powers 
of  c.  Tha  nxnrmatrix  f,  is  a  block- disgons 1-ms trix  of  the  form 


(1.3) 


f,(t»z) 


f*,(t,s)  0 


Tha  square-ma trices  f  ^  and  f*  ara  of  diawnsion  n_  and  n+  (where  n_  +  U+  •  n 
holds),  fj  has  only  strictly  stable  and  f *  has  only  strictly  unstable  eigenvalues  for 
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t  e  [0,1]  and  z  in  a  car tain  domain  D  of  ■*  (dafinad  aora  precisely  in  chapter  3). 
From  the  analysis  of  singular  perturbations  (see  c. f •  Hoppenstaedt  (1971),  O'Halley  (1975), 
Eckhaus  (1979),  Flaherty  and  O'Malley  (1980)  and  Ringhofer  (1981)  we  know  that  under  these 
assumptions  we  can  expect  a  smooth  solution  in  the  interior  of  [0,1]  while  at  the 
endpoints  boundary  layers  will  occur.  Therefore  the  usual  discretisation  methods  are  not 
applicable  unless  we  impose  conditions  on  the  mesh.  (See  c.f.  Eschar  and  Weiss  (1980)  for 
the  behaviour  of  collocation  methods  of  Gauss,  Radau  and  Lobatto-type. )  The  standard 
theory  for  discretisation  methods  for  general  grids  is  not  applicable  unless  the  maximal 
stepsize  is  smaller  than  e.  Flaherty  and  O'Malley  (1980)  avoid  this  difficulty  fay  solving 
only  the  reduced  problem  (which  is  independent  of  e)  instead  of  solving  the  full  problem 
(1.1):  As  e  tends  to  zero  the  solution  (y,z)  of  (1.1)  will  converge  to  a  function 
(y(t),z(t))  uniformly  in  compact  subsets  of  (0,1)  with  a  convergence  rate 
O(e).  obviously  (y,z)  satisfies 

(1.4)  0  -  f (t,y,z,0)  ,  s'  -  g(t,y,z,0)  . 

To  represent  (y,z)  as  the  solution  of  a  boundary-value  problem  (the  reduced  problem)  we 
need  m  boundary  conditions  for  the  equation  (1.4).  In  the  quaailinear  case  this  can  be 
done  by  a  "cancellation  law"  (see  Flaherty  and  O'Malley  (1980)).  In  compact  subsets  of 
(0,1)  this  gives  the  solution  (y,z)  up  to  the  global  discretization  error  we  make  by 
approximating  the  reduced  solution  and  a  term  of  order  0(e).  Inside  the  layer  this  gives 
us  no  approximation.  Using  some  information  about  the  analytical  structure  of  the  solution 
we  construct  a  method  which  is  designed  especially  for  problems  of  the  type  (1.1):  It 
shall  solve  (1.1)  directly.  So  we  need  not  compute  the  boundary  conditions  for  the  reduced 
problem  numerically.  The  amount  of  labour  (i.e.  of  gridpoints)  which  is  required  to  obtain 
a  uniform  approximation  of  the  solution  on  [0,1]  shall  be  "reasonable"  (that  means  not 
proportional  to  e  ').  If  we  use  a  general  (so  for  instance  a  uniform)  amah  with  maximal 
stepsize  H  >  >  e  the  approximation  shall  be  as  good  as  the  one  obtained  by  solving  the 
reduced  problem  only.  The  solution  will  consist  of  a  smooth  part  and  two  "layer-parts" 
exponentially  decaying  from  the  endpoints.  More  precisely  we  have 


y“(t,e)  -  y"(t>  ♦  u(|)  ♦  0(e) 
y*(t,e)  »  y+(t)  ♦  v(~^)  ♦  0(e) 

s(t,c)  -  s(t)  ♦  0(e) 

<y”  end  y+  ere  the  first  n_  end  lest  n+  components  of  y. )  u  end  v  decey 
exponentially  es  their  arguments  tend  to  infinity.  u<^)  behaves  roughly  like  the  solution 
C  of 

(1.6)  e  C'(t,e)  -  X  C(t,e),  c(0,e)  -  1,  Re  X  <  0  . 

A  discxet last ion- aathod  which  gives  us  e  good  approximation  of  y“  for  an  arbitrary  aesh 
should  also  be  able  to  deal  with  (1.6)  on  a  unifora  grid.  M  now  coapare  the  behaviour  of 
the  Box-scheae  and  the  implicit  Ruler  scheaa  at  the  problea  (1.6).  We  have 


(1.7) 


en,-cn) 

eLtl*L- 


1  (,«>  (1), 

2  l?i  Ci+l’'  * 


(2)  (2) 
cl+1  C1 


.(D 


(1 


2fc.i<#  W»,-i  .  ,.l  .(2) 

*  2?  '  Ci 


(2) 

1+1 


(1  -  -&)'1  ~  e 


lXh 
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(1.7)  shows  that  only  the  lnpllcit  Ruler  scheae  gives  us  a  good  approxiaatlon  of  the 
solution.  The  reason  for  that  is.  that  the  growth  function  of  the  Box-scheae 
(1+sXI-s)"1  tends  to  -1  as  s  tends  to  infinity  trtiereaa  the  growth  function  of  the 
iaplicit  Ruler  scheaa  (1-s)”1  tends  to  cero.  This  is  the  basic  idea  of  a  aethod 
developed  by  Xreiss  and  Richols  (1975)i  They  use  the  iaplicit  Ruler  scheae  for  y”.  Since 
the  eigenvalues  of  f*  are  strictly  unstable,  they  use  the  explicit  Ruler-schaae  for 
y*.  For  s  they  use  the  Box-scheae.  We  extend  this  idea  to  higher  order  methods  and  our 
approach  for  this  extension  is  polynomial  collocation.  For  a  positive  integer  q  we 
choose  q+1  reference  points  0  ■  u^  <..<  ■  1  on  the  interval  [0,1]  according  to 

(1.S)  u.  -  ^  ,  J  -  0(1 )q  . 

1  4 

For  a  given  aesh  T^  i “  (0  *  tQ  <..<  t^  ■  1,  H  ■  maxlt^-t^}  we  choose  the  collocation 
points  t^j  according  to 


(1.9) 


tii“ti  +  hiui  hi  *"  *1*1  *  V  3  “  o(,,q'  1  "  0(1  *,“1  • 

Note  that  ttq  »  ^t+1,0  holda'  0,1  *h*  subinterval  (ti#t1+1 )  we  now  use  for  y“  the 
collocation  points  tii«**,tiq‘  ,or  y+  w*  use  tiO'”,tl,q- r  ror  *  **  “** 

ti0,..,tiq.  For  q  -  1  this  gives  the  method  of  Kreiss  and  Nichols  except  that  we  use  the 
trapezoidal  rule  for  s. 

He  show  that  if  the  reduced  problem  has  an  isolated  solution,  the  operator  built  by 

P  and  the  boundary  conditions  is  stable  in  the  sense  of  Keller  (1975)  and  that  the 
ch 

stability  constant  is  independent  of  e  and  the  grid  TH.  So  if  we  use  a  mesh  selecting 

strategy  to  achieve  an  over  all  accuracy  6  we  only  have  to  control  the  local  discretize- 

tion  error.  He  develop  such  a  stepaize  algorithm  and  show  that  the  amount  of  gridpoints 

required  to  achieve  an  over-all-accuracy  S  is  essentially  independent  of  e  and 
_  2 

proportional  to  5  q.  This  is  comparable  with  the  Amount  of  gridpoints  our  method  would 
need  for  an  unperturbed  problem.  If  we  do  not  use  a  atepsize  algorithm  we  can  show  the 
following  result:  For  an  arbitrary  mesh  we  define  the  global  error  eh  by 
eh  "  (yh  -  y(th<£),  zh  -  z(th,e))  where  (y(th#e),  z(th,e))  is  the  solution  of  the 
continuous  problem  restricted  to  the  grid  V  Then  we  have 


(1.10) 


le^l  «  c[EHq+«<I+1-t-e+exp(-  ■—■J+exptX-^ — ) 

i-1  xh  -1  N-1  Kh  -1 

+  n  (1  ♦  ~~)  +  it  (1  ♦  —r>  ],  j  -  id )q»  i  -  0(1  )n-i 

r-0  €  r-i+1 


for  some  positive  constants  c,k,X  (here  H  denotes  the  maximal  atepsize).  So  away  from 
the  boundary  we  have  the  normal  order  of  convergence  0(Hq+1)  of  our  method  plus  a  term  of 
order  0(e).  Thus  the  approximation  is  as  good  as  if  we  had  solved  the  reduced  problem 
only. 

This  paper  is  organized  as  it  follows.  In  chapter  2  we  define  the  collocation  scheaw 
and  introduce  some  basic  notations.  In  chapter  3  we  present  some  analytical  results  of 
singular  perturbation  theory.  Por  the  proofs  we  only  refer  to  the  literature.  Furthermore 
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we  impose  the  main  assumptions  to  be  satisfied  by  the  continuous  problem.  In  chapter  4  we 
state  our  main  results  and  explain  them  In  a  sore  detailed  Banner.  In  chapter  6  we  prove 
thee.  In  chapter  5  we  demonstrate  our  method  on  a  model  problem  where  the  continuous 
solution  Is  known. 

2.  The  method. 

In  this  chapter  we  define  the  collocation  method  and  Introduce  some  basic  notations 
which  we  will  use  throughout  the  rest  of  the  paper. 

(2.1)  Definition!  We  call  a  set  of  q+1  points  0  <  uQ  <  Uj  <..<  <  uq<i  a 

reference  points  in  the  unit  Interval. 

(2.2)  Definition!  Given  N+1  gridpoints  0  “  tg  <..<  tj,  •  1  and  a  set  of  q+1  reference 
points  we  define  for  each  subinterval  Itjitj+)1  q+1  collocation  points  by 

(2.3)  tti  t±  +  h^j  hj  «-  ti+,  -  ti#  j  -  0(l)q,  1  -  0(1)11-1  . 

(Note  that?-  ti0  “  *  ti  holds.) 

(2.4)  Notation!  We  denote  the  neah  consisting  of  the  gridpoints  and  the  collocation  points 
and  its  maximal  stepslse  by 

Th  i-  {tlj#j  -  0(1 )q,  1  -  0(1 )N-1 },  H  -  max  hy  . 

(2.5)  Definition!  For  a  given  mash  TH  we  define  the  scalar  (vector-,  matrix-)  grid 
function  xh  as  a  sequence  of  numbers  (vectors,  matrices)  by 

*h  •-  {xijf  j  -  0(1)q,  1  -  0 ( 1 )N- 1 ,  xlq  -  xl+1  fl)  . 

(2.6)  Definition!  In  the  space  of  grid  functions  we  introduce  a  norm  by 

lxh«h  -  maxtlx^l^  j  -  0(1)q,  1  -  0(1)N-1) 

where  1*1  denotes  the  modulus  or  the  max-norm  in  rf1  or  the  matrlxnorm  introduced  by 
the  max  norm. 

(2.7)  Notation;  We  denote  the  first  n_  components  (rows)  and  the  last  n+  components 
(rows)  of  the  vector  (matrix)  x  with  x”  and  x+  where  n_  and  n+  are  the  dimensions 
of  f"  and  f*  in  (1.3). 
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>nt  To  apply  a  (polynomial)  collocation  method  on*  ganarally  seeks  a 


(2.8)  Deflnltlo 

polynomial  spllna  function,  which  aatiafiaa  tha  given  differential  aquation  at  certain 
points  (tha  collocation  points)  as  well  as  tha  boundary  conditions.  In  order  to 
approximate  tha  solution  <y<"  (y”,y+)T) , z)  of  (1.1)  we  construct  a  vector-spline  function 
<Py»Pz>#  Py  "  (Py«Py)T  which  satiafiest 
a)  p  ,p  e  c,.  ,, 

*y  *z  10,1) 

8)  Py  and  p^  are  polynooLials  of  degree  qt  pz  is  a  polynomial  of  degree  q+1  in  each 
aubinterval  (ti(ti+1) . 

Y)  To  generalise  the  implicit  and  explicit  Euler  method  we  request  that  p~  satisfies  the 
differential  equation  at  all  collocation  points  except  the  first  one  in  each 

subinterval  and  p*  satisfies  the  differential  equations  at  all  collocation  points  except 
the  last  one  in  each  aubinterval. 


(2.9) 


e(p  )  -  f  <t,p  ,p  ,«)  for  t  -  t  .,  J  -  1 ( 1 )q,  1  -  0(i)H-1 

y  y  * 

e(p*)'  -  f+(t,p  ,p  ,e)  for  t  -  t. .,  J  -  0(l)q-l,  i  -  0(1  )N-1  . 

y  y  *  *j 


Por  z  we  take  all  q*1  collocation  points  except  for  2  subintervals  .  For  reasons  which 
will  be  explained  in  chapter  6  it  is  necessary  to  modify  the  method  used  for  s  in  order 
to  obtain  a  higher  order  method.  Given  a  mesh  TR  we  assume,  that  I,  steps  h^  on  the 
left  hand  side  and  I2  steps  on  the  right  hand  side  ere  of  order  0(  c)  (where  we  also 
allow  If  and  I2  ho  be  zero): 

ht  -  o(*>  i  -  o(i)i1  -  i,  i  -  i2om-i 


h.  >>  e 

1 


U-I^l 


»  e 


For  the  subinterval  [t-  ,tT  +1]  we  take  the  collocation  points  tT  ,,..,tT  _  (as  for 

*1  1  ‘v? 


Py).  For  'Vle-vVlJ 


we  take 


tH-I.-1,0***'tll-I~-1'  q  ”  1  *°  h*T*‘ 


y  '-.ii 


‘Wij 


iO 


(shV 


ij 


-  ? 

3-1 


jk  lk 


(ShVtj 


*1 

j-o 


jk  lk 


<Sh*h)13 


f'Wij  if 

l  KVij  tf  1  "  ,"I2-1 


%  0 

wjk  *lk 
k«0  ■* 


With  this  notation  tha  discrete  problea  can  be  written  as 
(2.7)  Peh(W  '  ° 


P 


Ch'W 


6  Vh 


6  Dh  yh 


D  Z 
h  h 


sh  f  <th'yh'Bh'e> 

•  sh  f*(vvve> 
<  «<vwc>  • 


(2.8)  Rsurki  To  this  point  we  did  not  restrict  ourselves  to  a  certain  set  of  reference 
points.  1  natural  choice  would  be  Lobato- points  in  order  to  achieve  the  highest  possible 
order  of  accuracy  for  z.  However,  for  the  rest  of  this  paper  we  restrict  ourselves  to 
equidistant  reference  points. 


J 


3.  Analytical  results  and  main  assumptions 

Quasilinear  singularly  perturbed  b.v.p.'s  of  the  for*  (1.1)  have  been  studied  by 
several  authors  (see  c.f.  O'Halley  (1974),  Flaherty  and  O'Malley  (1980),  Howes  (1980), 
Ringhofer  (1981)  ).  In  this  section  we  present  their  results  in  a  manner  which  is 
appropriate  for  our  purposes.  One  can  find  the  proofs  either  directly  in  the  above 
mentioned  papers  or  can  prove  the  results  of  this  section  analogously  using  the  special 
structure  of  problem  (1.1). 

The  solution  (y,z)  of  (1.1)  has  a  uniformly  valid  asymptotic  expansion  of  the  form 

«0 

(3.1)  y(t,e)  ~  l  <?.( t)  +  y.(^)  +  v .(■— Oe3 

j-0  3  36  36 


(3.2) 


m  c» 

*(t,e)  ~  l  z  ( t > £3  *  l  (o  (|)  + 

j-0  3  j»1  3  3 


where  are  exponentially  decaying  as  their  arguments  tend  to  infinity.  (We 

call  a  function  4»< T )  exponentially  decaying  If  it  satisfies  M>(t)I  <  c1  exp(-c2r)  for 
positive  constants  c-j  and  c2.)  To  derive  this  expansion  it  is  necessary  to  define  the 
reduced  problem.  This  is  done  in  Flaherty  and  O'Malley  (1980)  for  the  case  of  separated 
boundary  conditions.  In  our  case  we  proceed  analogously. 

(?0'V'M0  “nd  vo  satisfy 

(3.3)  0  “  f (t,y0,zQ,0 ),  z'  ■  g(t,y0,zQ,0) 


(3.4) 


f?  VT)  "  fi(0'*o(on*‘,o,T)'  I?  vo(T)  “  -*1«i.?»D)v0(t) 


(3.5) 


I  ?o(0)  +  V°A 

b  (z  <0),z  (1))  I  +  b  (0,0,z„(0),z  (1),0)  -  0  . 

\y0(1)tv0(0>y 


Since  vQ  and  u0  must  decay  exponentially  Uq(0)  “  0,  v”(0)  “  0  must  hold, 
up  b^ ,  into 


(3.6) 


b1(W 


b11<‘0'*1> 


t  *12<VV 


We  split 
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(Htr<  y**  and  y*  ara  dafinad  as  in  So,  bf  j  and  bj2  (n+*^  *  n 

2ib 

matrices.  We  assume  that  b11  has  maximal  rank  in  a  domain  of  R  containing 
(*0(0),r0(1))  and  that  there  exists  a  regular  (n+m)x(n+m)  matrix  WtZg.z^ )  much  that 


(3.7) 


»<V‘l  *b11 **0'*1  * 


holds  where  V  is  a  regular  nxn  matrix.  In  that  ease  (3.5)  is  equivalent  to 

’~f  \ 

y0,0> 


(3.6)  W  (r0(0),*0(1))(b12(*0(0),x0(1)) 


y_(U 


+  b2(0,0,E0(0),x0(1),0>)  -  0 


(3.9) 


rV0> 

lV+(°) 


-  -v”1  (z(0 ) ,z( 1 )  )w”(z(0 ),z(  1))b(y(0),y(1),s(0),z(1),0) 


(Here  w"  and  W+  denote  the  first  n_  and  the  last  n+  rows  of  W. )  (3.8)  together 

with  (3.3)  gives  the  reduced  problem.  If  the  reduced  problem  has  an  Isolated  solution  we 
can  show  the  uniform  validity  of  the  asymptotic  expansion  (3.1),  (3.2).  This  leads  to  the 
following  assumptions :  There  exists  a  domain  D  of  Rm  so  that 
(3.10)  HI i  f  ^ ( t , x )  is  a  block  diagonal  matrix.  So 


f^t.z) 


f*(t,z) ; 


holds.  There  exist  positive  constants  X_  and  A+  so  that  all  eigenvalues  1  of  the  n_ 
dimensional  block  f”(t,z)  satisfy  Re  X  <  - A  and  all  eigenvalues  of  the  n+  dimensional 
block  f*(t,z)  satisfy  Re  A  >  A+  for  all  t  e  10,1)  and  s  e  D. 

(3.11)  H2«  There  exist  matrices  WtZg.z^  and  VfZg.Zj)  as  defined  in  (3.7)  for  all 
C  De 

{3.12)  H3i  The  reduced  problem  (3.3),  (3.8)  has  an  isolated  solution  (y0(t),z0(t) ) . 

Zg(t)  lies  within  a  compact  subset  of  D  for  all  t  6  (0,1). 

(3.13)  H4»  f,g,  and  b  are  as  often  differentiable  as  necessary  with  respect  to  all  of 
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•>  •*  • 


"i* 


their  arguments. 


L 


Under  thee*  assumptions  the  asymptotic  expansion  (3.1),  (3.2)  is  uniformly  valid. 
(3.14)  notation i  For  the  further  we  denote  with  (Wb)(y,r)  •  0  the  boundary  condition. 
H(s(0,c),  s(  1,  e)  )b(y(0,  e),y(  1,  e),*(0,  e),s(  1,  e) )  -  0  . 


4.  Main  result* 

if*  f ^r>t  show,  that  under  the  assumptions  (3.10)-(3.13)  there  exists  also  a  solution 
of  the  discrete  problem 

«•’>  reh(yh'V  “  °*  be(yh'V  “  0 

and  this  solution  is  stable  in  the  sense  of  Xaller  (1975).  Me  start  with  constructing  a 
uniform  o(c+H)  approximation  of  (y^s^i  He  define  the  discrete  reduced  problem  by 

(4.2) 


rh<yh'V  “  °*  1>R<yh'*h>  “  ° 


/  rtti3'yij'*i3'0> 

(P",  .  j  -  1(1  )q,  i  -  0(1)N-1 

WjVV^WV0” 

bR(yh'V  "  W2(*O0'*M0>b<y00'yM0'*OO*‘HO'0)  * 

Because  of  hypothesis  (3.12)#  (4*2)  has  an  isolated  solution  (•••  Beller 

*  0  0  . 

(1975)).  We  now  define  yh^h'*h  hjf 

^  ■  VVVVV*  mhyh  ■  VVVVV 


, 

,5“  •] 


(4.3) 


*+  ]  “  “V^*00'*H0^  W1<*00'*K0,b<y00'yHO,*0O'*H0'°) 


rN0 


0 


0 


W  ‘h  “  lh  • 


(y^,*^)  is  a  uniform  0(e)  approximation  of  the  solution  ( y^ » *1, )  of  (1.18).  To  show 

the  existence  and  the  stability  of  (yh,*h>  we  need  that  the  linearisation  of  the  operator 

0  0 

built  by  Feh  together  with  the  boundary  conditions  (3.14)  at  (y^,*^)  *n  inverse 

bounded  uniformly  in  e  and  H.  This  gives  a  linear  difference  operator  (L^iB)  of  the 
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i 


form 


(4.3) 


(4.4) 


Dhyh  "  8hl(Ah  1  yh  +  ,Ah  )  *hA 


WVV  ‘I  e  Vh  ‘  8ht(*h1,+yh  + 

„  0 .  21  .  .23  , 

Vh  •  V\  yh  +  Ah  V 


/It  -  -f  T  t2  +  ~  T 

8  ,<y  OO'W  *  8  {y0O'yMO'*OO'*H0) 


8(W 


\8  {y00'yN0',00'*H0) 


Sine*  y8  will  exhibit  a  rapid  tranaition  at  tha  andpointa  tha  aatricaa  A^ 
of  a  aaooth  part  and  a  part  exponentially  decaying  away  froai  tha  boundary  i.a. 
satisfy 

.ki  _  r*i  *  :it<  ,n  c*i,  < .  w  ...  ,  .  - 

\  Ah  \  '  ’Vh  'h  -  V*'1  1,a#  Ah  ° 


(4.5) 


*21 
Ik  I 

ij 


1-1 

Oh 

H-1 

oh 

n  (i 

♦  ♦ 

n 

<1  ♦ 

r«0 

€ 

r-i+1 

1-1 

C2  " 
r-0 

*hr  -1 
(1  ♦  -f) 

HA 

12  «■  M_1 
15>  '  5  " 

3  r-l+1 

Yh 


k,-i 


j  -  0(1)q,  i  -  0(1)H-1 


will  consiat 
they  will 


Because  B  is  the  linearisation  of  the  boundary  conditions  (3.14)  B  will  satisfy 
detlB11)  t  0. 

For  the  problaa 

(4.6)  Leh(VV  “  V  8(W  "  8 

(where  ♦.  is  a  vector  grid  function  of  appropriate  dinension  and  0  e  Rn+n)  we  define 
n 

the  reduced  problem  by 

(4.7)  -  V  “X'V  ‘  ®2 


12 


(4.8) 


^‘vV  ■  { 


*  K2)\] 


■sh[(xh1,+yh *  K2)\] 


„  021  .22  , 

Vh  •  V\  yh  ■  Ah  V 


(4.9) 


B  ‘VV  "  B‘(y00'yH0**00'W 


(B2  denote*  the  last  a  component*  of  B).  For  the  ao  defined  linear  operator  we  have 
the  following  stability  result 

(4.10)  Theorem:  bet  L  .  ,B.LB.BR  be  defined  as  in  (4.3).  (4.4),  (4.8),  (4.9).  If  there 
Cn  n 


.R 


,(yh'*h),h  s  ci,,Lh(VIh),h  +  ,BR(VIh),) 


There 


exists  a  constant  c1#  so  that  satisfies 

(4.11) 

then  the  problea 

1 4 * 1 2 >  Leh(yh'V  "  V  B(W  -  6 

has  a  solution  ( yh , z^ )  for  all  vector-grid  functions  4^  and  all  vectors 
exists  a  constant  c2  independent  of  e  ao  that  (yh,*h>  satisfies: 

(4.13)  "VV,hiC2"U+'B')  * 

Since  the  continuous  reduced  problea  has  an  isolated  solution  the  linearization  of 

(4.14)  yeh<W'  <WbHyh''h> 

satisfies  the  hypotheses  of  theorem  (4.10).  (Here  (Wb)(yh,zh)  is  defined  analogously  to 
^3.14).)  Using  theorem  (4.10)  we  can  show  the  existence  of  a  solution  of  (4.1).  Moreover 
the  operator  defined  by  F£h  and  the  boundary  conditions  wb  is  stable  in  the  sense  of 
Keller  (1975)  in  a  neighbourhood  of  (y^.*!,)  of  (4.1).  The  neighbourhood  and  the 
stability  constant  are  independent  of  e.  So  we  have 

(4.15)  Theorem:  If  the  continuous  problea  (1.1)  satisfies  the  hypotheses  in  chapter  2  then 
there  exist  some  positive  constants  cq>hq  *°  that  (4.1)  is  solvable  for  all  e  <  cQ 

and  H  <  Hq.  The  solution  f of  (4.1)  satisfies 

(4.16) 


l(yh,sh>  -  (y°,s8)lh  »  °3 * 


0  0. 


for  soae  positive  constant  c3  independent  of  e.  Here  (yh,«h>  are  defined  as  in 
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(4.3).  Tha  atablllty-lnaquallty 


(4.17) 


,<yh'*h)_(yh'*h,'h  i  C4,F€h‘V*h,"F«h(yh'*h,,h  + 


♦  l(Wb)(y^,a^)  -  (MbXy^,*^)  I 


la  aatiafiad  for  all  i  "  l.a  aatiafying  KyJJ#*^)  -  “a  *11 

«  <  V  *  *  V 

Sine*  (y®,«^)  In  (4.3)  la  an  O(e)  approximation  of  (y^ty,)  tha  global  arror 

<<y^.«^>“<y*t.  ,*),  alt.  ,*)))  la  up  to  tana  of  ordar  O(c)  glvan  by  tha  diffarance 
n  n  n  n 

batwaan  (y°,a°)  and  tha  0(1)  tarma  In  tha  aayaptotic  axpanalon  (2.1),  (2.2)  of 
(y(t,*),a(t,t) ): 

(4.18)  Thao  ram:  Lot  (y0,*0)  danota  tha  0(1)  tarrn  In  tha  aaymptotlc  axpanalon  of 
(y(t,e),*(t,e) ).  Than  thara  axlat  poaitlva  conatanta  c  and  tc  ao  that 
,(y0(tiJ,*).a0(tij.O)  -  (yJj.aJj)!  <  cW*^1  ♦ 

1-1  xh  »-1  *«h  T.t 

R  (1  +  —*)  +  H  (1  +  -~)  ♦  axp( - * 

r-0  *«i+1 

tlL^ 

•xp(y+  -«*— )) 

holda  for  all  H  <  H„  and  «  <  cQ  whara  T_  and  r+  ara  tha  (In  nodulua)  aaallaat  raal 
parta  of  fj(0,»(0))  and  f*(1,s(1)). 

From  Thaoraai  (4.18)  followa  lmmadlatalyi  ^ 

(4.20)  Thaorami  Lat  (y(t,c),a(t,c))  and  (yh,a^>  ba  tha  aolutlona  of  (1.1)  and  (3.1). 
Than  thara  axlat  poaitlva  conatanta  cfi,  y  and  <  ao  that 

l(v(t.  ),*<t.  .))  -  (v.  . )  I  <  cXH*1*1  ♦  e  ♦ 


■*r-’ 


•ly(tij),*(tij ))  -  (yAJ 

1-1  ah  , 

M-1 

n  (i  *  t*)  ♦ 

s 

r-0  * 

r-ial 

t  -1 

axpd.^-)) 

holds  for  all  c  <  c0  and  all  grids  with  H  <  Hq. 


Away  tram  the  boundary  wo  will  hawa  an  approximation  of  ordar  Ott1,  +e)  for  an 
arbitrary  (ao  for  instanea  alao  unifora)  nosh.  So  in  this  caaa  wa  hawa  tha  ordar  of 
conwergance  which  wa  would  obtain  by  solving  tha  raducad  problaa  plus  a  tara  of  ordar 
0(c).  In  that  sansa  our  aathod  is  equiwalent  to  solwing  only  tha  raducad  problaa  axcapt 
that  wa  naad  not  coaputa  tha  boundary  conditions  of  tha  raducad  problaa.  Wa  now  oonsidar 
tha  casa  wharo  wa  want  a  battor  approximation  than  tha  ona  wa  obtain  by  using  an  arbitrary 
sash.  This  is  of  saas  interest  if  either  c  is  rather  big  (let's  say  c  -  10~2)  or  if  we 
are  interested  in  tha  solution  ins ids  tha  layers. 

If  wa  can  show  now  that  (y(th,c),s(th> e) )  (tha  restriction  of  tha  exact  solution 
(y(t,e),s(t,c))  on  tha  aash  Tg)  lias  within  tha  ball  Kp(y4,*4).  thsoraa>  tails 

us  that  we  only  have  to  control  tha  local  discretisation  error  to  obtain  a  unifora  approxi¬ 
mation  of  tha  solution.  In  other  words  if  wa  choose  tha  wash  TR  ao  that 

wA,y(v«»'  *(ve,,,h  < 4 

holds  wa  hawa  l(yg#s^)  ”  <y(hjl»*)»«(th,c) ) I  •  0(8). 

(4.22)  lawai  Lot  tha  local  discretisation  error  (**,|£)  be  defined  by 


(4.22) 


“S'#  •FC.hWV#>'-,Vrt> 


where  (y(th,c),s(th,c))  denotes  the  restriction  of  tha  solution  of  (1.1)  to  tha  grid 

Tr.  Than  there  exists  a  posltiwa  constant  c7  such  that 

,  h  1  t 

*  “ij'4ij  1 1  •  c7<*hi  4  hi  4  (e,q  •XP<"  -^r> 


(4.23) 


*  ■j4-'*’  — - 


h  g  t  -1  h  q+1 

4  <?>  ir> 4 


Tq^lTc” 


holds  for  J  -  0(1  )q,  i  -  0(1)11-1  and 


ti  '“i'1!*!** 


So  if  wa  want  to  guarantaa  that  Kfj^, #  *h  "  holds  for  a  certain  desired 

accuracy  (  we  hawa  to  choose  tha  mesh  T_  according  to 


-  -  -  -2-  -  X  t  -L~ 

h^  <  ain{4q  E  q,  4q*\  4q  €  exp(-^“),  tf***  E 


(4.24) 


X  t. 


-X.  (t  -1 ) 


A  t.  —  —  A  X  U  .  -  •  f  "TT  "^+^4*^  ^ 

•xp  7^777  •  5  e  “V'—C - »'  **  e  exp<  (^TIT-'1 


If  we  choose  according  to  (3*24)  we  can  show  that  My^# )*(y(t^# t) ,z(t^# €)  i  * 

0  0 

OU+6)  holds.  Therefore  (y,s)  lies  within  the  ball  with  radius  p  and  center  <y h*zh> 

where  F  is  stable  (see  Theorem  (4.15))  if  only  e  and  6  are  sufficiently  small, 
eh 

Applying  Theorem  (4.15)  gives 

(4.2S)  Theorem:  If  TH  is  chosen  according  to  (4.24)  then  there  exist  positive  constant. 


c8,  e0  and  6^  so  that 

(4.26) 


,<yh'“h>'{ytth'E)'Z<th'£,),h  ‘  C66 


holds  for  any  6  <  6Q  and  E  <  Eg. 

To  construct  the  mesh  TH  we  proceed  as  it  follows: 


x-v 


Step  1:  Starting  with  h0  we  choose  hA  according  to  <  min  { edf5  exp( 1  ~~)  , 


e«q+1  exn  — as  long  as  h.  <  6q+1  holds.  (For  a  reasonable  4  >  cq+1 
(q+ 1 )€  1 

i  i  2 

6q+1  <  6q  eq  holds.) 

Step  2:  Starting  with  h^.,  we  choose  hi  according 

.  .a  +  i+1  .0+1  +  i+1  , 

hi  <  min{c6  exp  - - - ,  e4’  exp  - TTTTT  J 


he 


(q+1 ) E 


as  long  as  hA  <  4q+1  holds. 


.q+1 

Step  3:  Por  the  rest  of  the  interval  choose  <  0’  . 

Therefore  it  Is  necessary  to  know  the  constant  X_  and  X+  which  are  the  smallest 
moduli  of  the  real  parts  of  f”(0,r(0))  and  f*(1,*(1)).  These  can  either  be  computed  by 
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solving  ths  problem  on  *  uniform  stash  first  or  eatiamtes  tor  than  can  ba  obtained,  if  ona 
uses  a  continuation  method. 

An  * - dlately  arising  question  isi  How  big  is  the  additional  mount  of  labour 

insids  ths  layers?  That  means t  How  many  grid-points  do  we  need  in  steps  1  and  2?  What 
one  would  like  to  have  is  an  estimate  for  the  number  of  grid-points,  which  is  independent 
of  c.  Unfortunately  the  number  of  grid-points  required  to  achieve  a  certain  accuracy  S 
tends  to  infinity  as  c  goes  to  zero,  but  the  growth  is  so  slow  that  for  reasonable  ranges 
of  c  and  4  it  appears  to  be  independent  of  e.  More  precisely  we  have 

(4.27)  Lesmmt  Let  N(£,X,y,c)  denote  the  number  of  grid-points  produced  by  step  1  (or  2) 
(so  X  either  equals  X_  or  X+) .  Than  we  have 

(4.28)  |N(£,X,q,C^)  -  N(5,x,q,t^)|  <3 

for  6,X,q  and  and  c2  satisfying 


Furthermore 


.  -18  ,  ,  ..-2 
10  *  Ei  *  10 


i  -  1,2  . 


-  1 

•  (4.29)  n(£,X,q,c)  «  const*  6  q  *n(ln  (^)) 

holds. 

(4.30)  Remark «  The  estimate  (4.29)  is  not  sharp,  but  together  with  (4.28)  it  shows,  that 
for  reasonable  choices  of  t  and  £  the  amount  of  additionally  needed  grid-points  is  of 
the  same  order  or  magnitude  as  a  q-th  order  collocation  method  would  need  for  on 
unperturbed  problem. 


5.  numerical  Example 

We  now  illustrate  our  results  with  a  third  order  quasilinear  model  problem  where  we 
know  the  exact  solution.  For  the  computations  a  modified  version  of  the  packages  LOMTO 
and  SOLVXBLOCX  (see  de  Boor-Weiss)  (1976))  has  been  used,  we  consider  the  system 
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*  r 


y1(1<c)'fs(1,c)-2«  -  c 


Hara  tha  aatrix  la  of  tha  for* 

(5.3)  f,(t..)  -  r*  ®) 


■oti  that  tha  arror-conatanta  ara  dacraaaing  for  incraaaing  q.  (00  for  for 
q  “  4  the  error  constant  is  0(10-2)).  To  achieve  tha  accuracy  tha  following  Mounts 
aeshpoints  have  bsan  required) 


Required  eeshpointsi 


i\« 


'  \ 

1 

2 

3 

4 

25 

17 

25 

33 

41 

13 

26 

43 

58 

74 

63*- 1 

44 

77 

109 

140 

311-1 

78 

144 

208 

272 

•o  Cor  a  Clued  q  the  auaber  oC  required  eeshpolnts  grows  like  «q.  Varying  t  Crow 
10  to  10  the  awount  oC  grid-points  Cor  Cixed  q  end  4  varyiad  only  by  ±1. 

How  we  solve  (4.1),  (4.2)  using  s  unlCore  wash  (hA  ■  H,  1  ■  0(1M-1).  Because  oC  the 
theory  in  the  preceding  chapter  we  expect  an  error  oC  order  (l**1**)  away  Crow  the 
boundary  layers.  Xn  the  next  table  the  norn  o C  the  error  at  t  -  0.5  is  listed  Cor 


m- 


.1 

.51-1 

.251-1 

.11-2 

.581-3 

.551-3 

.511-3 

.11-4 

.531-5 

.  61-5 

.  61-5 

.11-6 

.  51-6 

.  41-7 

.571-7 

So  we  aee  that  the  error  la  of  order  c  Independent  of  H  except  for  e  -  .  IE-6  end 
H  ■  .1  where  cHq+1  >  e  holds .  <c  is  the  error  constant  of  our  method  for  q  “  3. ) 


6.  Proofsi 

To  prove  Theorem  (4.10)  we  first  need  soae  auxiliary  results.  He  start  with  scalar 
constant  coefficient  probleas  and  estiaate  the  growth  function  of  the  aethod  we  use  for  the 
coaponent  y”  in  (1.18)< 

6.1  La— ai  Given  the  complex  numbers  Xg  and  s  and  the  vector  e  ”  (e1 , .  •  >*q)  the 
solution  x  “  (x1,...,x^)  of  the  systaai 


(6.2) 


Kj  '  *0  +  *  Vk  +  '  3  ■ 

satisfies 

Xj  -  Yj(«)x0  ♦  6j(s)e  . 


are  vectors  of  dimension  q.  y^  and  the  components  of  B^  are  rational  functions 

of  s.  For  q  ■  1(1)10  there  exist  angles  8  satisfying  0  <  0  <  x  and  positive 

q  q  * 

constants  c_  such  that 

(6.3)  |y,(z)|,  IB, (z)  I,  Izy.  (z) I .  IzB.  (z)  I  <c 

j  J  j  3  q 

holds  for  s  in  8  ■  (s|Rez  <  0,  I  la  z(Rs  *)  “ 1 1  <  tan  8  }  and  j  “  KDq.y  satisfies 
q  q  q 

It  I  <  1  in  8  and  y  (0)  -  1. 
q  q  q 

Proof i  (5.2)  can  be  written  in  the  fora 

x  -  "?*  Xg  ♦  zQ”x  ♦  e 

— ♦  f  •  a 

where  1  denotes  the  vector  (1,..,1)  and  0  the  q  x  q  matrix  (w^.i, J"1,q).  if 
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M  solve  this  syaten  by  eraser's  rule  and  look  at  tha  Involved  determinants  as  functions 
of  a  wo  can  writs  ths  solution  (x,,..,xq)T  in  ths  for*  (6.2).  Ths  and  0j  turn 

out  to  bs  rational  functions  whars  ths  degree  of  ths  denoMrator-polyncaial  is  q  -  1  and 
ths  degree  of  ths  dsnoailnator  polynonial  is  q.  Thorsfors  thsy  bshsvs  liks  s  1  at 

s  “ 


Ths  growth  functions  Y^  bars  boon  analysed  nunarically  by  Oabsrhubsr  (1979)  and  ars 


lass  than  1  in  Modulus  within  ths  sector  Sq. 


(6.15) 


t  Let  all  eigenvalues  X  of  the  nkrt  ns  trices  lie  within  a  oonpact 


subset  of  tha  interior  of  8q  and  satisfy  ReX<-X<0  for  a  positive  constant  X.  Let 


ID  <  c11  hold  for  a  positive  constant  c1f.  Then  ths  solution  yh  of 


(6.16) 


“h^^Vh’-'h  '  W  “tUfis* 

,yh'h  <  *n'M  +  ,fhV 


for  sans  positive  constant  c12  independent  of  t.  (Here  yh  and  fh  are  vector  grid 
functions  of  dinension  n.) 

Proof!  We  rewrite  (6.16)  in 

“Vh’iJ  -  »i9«Vh>lJ  *  V8h((°hVyh»iJ  +  fiJ  3  -  1  ’  • 

Since  Djjkjj  is  bounded  we  get,  by  using  lean  (6.1) 

h.  q  h,  h. 


'1+1 


.0  -  xi0)yi0  *  j,  <V^  ^  V  +  0(Vyh  V 


i 

Here  the  rational  functions  Y  and  0  are  applied  to  the  nonelngular  Matrix  —  A 

q  q  c  ID 


the  usual  wanner  i.e. 


-/t  y«^ 


where  f  is  a  closed  curve  in  the  coaplax  donain  and  all  eigenvalues  of  Ai0 
r.  (See  c.f.  Dunford  Schwarts  (1957).) 


lie  within 


(6.17) 


'  V 


M  N  rr  h  -1  h. 

Iy.„l  <  lal  +  const (  J  S  (1  -  — - — )  —■  * 

10  j-0  r-j  *  * 

<6,,7)  i-1  i-1  icy_»»  -1 

♦  l  n  (I---2)  )(lf  lh  ♦  0(Hly  I  )) 
j^+1  r-j+1 

h  -1  -itth  w 

by  estimating  <1  -  *)  <  exp(  ~~~)  with  SOM  suitable  constant  w  ws  gat 


MM  *y  h  *1  h.  M  -KYh, 


n  n  n  |  it  t  11  ,  n  n  |»  h  ,  .  tM  m  ■( 

(6.18)  I  I  (1  -  —p-)  -1  <  l  sxp(— P  wlt-tlp  <  i  /."sxpf-  ~T(— 7 - Ida  . 

j-0  t-j  *  6  j-0  *  m  j  e  e  o  e 

Sines  tha  intagral  is  boundad  and  tha  other  sum  in  (6.17)  is  geoMetric  wa  gat  the  dasirad 

estimate 

(6.19)  Laaami  Let  A^  be  a  n*n  Matrix  and  ?h  and  be  n-dinansional  vector  grid 

functions.  Lat  A^  satisfy  tha  conditions  for  tha  Matrices  A^  in  baana  (6.15).  Wa 


aasuae  that 


,DhVh  <  C13*  ,Dhfh'h  *  c13 


(6.20)  .  i  xh  -1 

lf..l  *  c14  "  (1  +  ~T)  1  “  HD**  1  “ 

3  r-0 

holds  for  sane  positive  constants  c13,  c14  and  x.  Than  tha  solution  yh  of 

(6.21)  “  °'  y00  *  «  **tl*flM 

.  .  1  oh  -1 

y.  -  -(A-1f>  ay*  0(e),  ly  I  <  c  H  (1  ♦  p> 

(6.22)  h  h  h  13  15  r-0  6 

j  -  1  (1  )q,  i  -  0I1W-1 


(6.22) 


for  som  positive  constants  c^g  and  a  whare  a  <  mln(x,c12)  holds.  (Bara  c12  ia 
defined  as  in  Lsmms  (6.15). 

Proof i  We  define  y^  and  u^  by 

(6.23)  CDhyh  -  VWV*  1  "  1<1,q'  1  “  y00  "  *  +  *00*00  ' 


(6.24) 


^  -  yh  a  (A  f)h  -  yh,  than  satisfies 


<6-25>  %%  ■  vw  ♦  ■vvv 


i .r*»  ;t  tr'^-Kr  -  .v> 


jw'-'V  >*  i  ■■■  '  f  ‘  ’  v  .  -t&i 


Because  of  Learns  (6.15)  v*  obtain  that  lu^  -  Ole)  holds.  Ha  define  »h  by 


(6.26) 


J*  °  *  1  “  . . 00  ”  y00  ' 


for  km  positive  o  leas  than  min(x,c6  ).  vh  aatisfiea 

1  °h 

(6.27)  eoAJvh  -  ovi0  +  pij  “  *  11  +  “7^  3  ”  KDq#  l  *  0,1,..  . 

Because  of  (6.20)  4  Ct4  holdB*  80  Th  lB  bound*d  and  therefore  (6.22)  holds. 

(6.28)  Lesmai  Let  the  operator  associated  with  the  problem 

<6*29)  Dh*h  “  Bl*oo'W  ■  n 

have  a  uniformly  bounded  inverse  i.e.  let  1*^  lh  <  c^dnl  ♦  1^)  hold  for  some 

constant  c^.  (Here  s^  and  denote  vector  grid  functions  of  dimension  s.  H,  n 

and  Cj,  are  of  appropriate  dimensions . )  Let  the  a- dimensional  vectors  fh  satisfy 


(6.30) 


lftjl  <  c1?(leljl  ♦  J  *  0(1 )q,  i  -  0( 1 )N 

i  c  h  -1  H-1  c.  h  -1 

e.  -  n  (1  ♦  — ,  j  -  1(1)q,  e..  -  n  (i  ♦  -1S-Z) 

11  r-0  *  r-i  ' 

j  -  0(1  )<j-l,  1  -  0(1)H-1 


for  soan  positive  constants  ci;  and  c1(.  Then  the  solution  of 

(6.31 )  ®h*h  -  8h<Ch*h  *  V  *  °'  "‘•oo'So’  '  Cfl 


(6.31) 

(6.32) 


'Vh  «  c19*^  +  c20> 


for  some  positive  c19  independent  of  c. 

(6.36)  leii  Let  L(h  be  defined  as  in  (3.3).  Let  the  coefficient  matrices  A^ 
satisfy  (3.S ) .  Assume  that  the  operator  associated  with  the  reduced  problem  (3.7)  has  a 
uniformly  bounded  inverae,  i.e.  that  there  exists  a  constant  cM  so  that 


(6.37) 


“W'h  «  e22<n‘h‘V*h,'h  +  •»R<yh^,,) 


holds  for  all  grids  TR.  Then  there  exists  a  solution  (yh,s^)  of  the  problem 

(6-38)  fcth(W  ■  V  y00  “  V  <0  “  V  bR(VV  “  #2 

T 

for  all  ( n+m)-deminslonal  vectors  B  -  (BQ,B^,^)  and  all  grid  functions  4^.  (y^s^) 

satisfies 


m 


(6.39) 


,(VV  *  C23l  *  'W 

(or  a  positive  constant  c23  independent  of  c. 

Proof i  Let  tha  n  *  n  matrices  Th  ba  dafinad  by 

•vv"  -  s;;«AhV(Thr  +  ,xh2,'}  -  °'  <too>"  -  ° 

(6.40) 

*VV+  -  S>XV(V+  *  <^l2>+,  "  °'  (tho)+  “  0  • 

(Kara  t"(  (A^'  )",(A^ )”  >  *nd  Th<<Ah1,+'  {A12)+>  tha  first  n_  and  last  n*  rows 

of  T.  (a!)4  ,A^)  ) .  Bacauaa  of  Lemma  (6.19)  aatlaflaa  T.  »  -(a"). 'll'*).  ♦  T  ♦  0(c) 

h  1 1  1 2  n  h  h  hi 

and 

1  ich  -1  .  N-1  ch  -1 

IT"  I  <  n  (1  ♦  — *■)  .  j  -  Id  )q,  IT  I  <  H  (1  +  -~) 

(6.41)  13  r-0  3  r-1 


j  -  0(1)q-1,  i  -  0(1  )K-1  . 

Ha  daflna  uh  -  yh-Thrh  and  obtain 

eoijuh  ■  -IXV *  *  hAiUKl')t<  *  "12>!vh!  - 

<6-«>  Vh  •  «\ ♦  K2  -  ♦  M'w  *  < 

uoo  ’  *o*  v  ■  v  ■  02  •  B"(wv  • 

So  for  c  and  H  sufficiently  small  (6.42)  can  be  solved  by  using  contraction  and  (6.39) 
holds. 


With  the  aid  of  Leaau  (6.36)  we  now  construct  the  general  solution  of 

L  (y.  ,z.  )  “  4  i.e.  a  solution  which  depends  on  ( n+m)  free  parameters.  He  start  with 
cn  h  h  h 

the  homogenous  problem i 

(6.43)  Lemma i  Let  tha  assumptions  of  Lemma  (6.36)  hold.  Then  there  exist  solutions 

(yj,«^),  ^h'*h*  of  th*  Probl,l“ 


(6.44) 


for  any  n 


Lch(y^lh)  ■  °'  ,yoo’’  ■  T1"' 
Sh'vs’  ■  °» lySor  ■  °' 

dimensional  parameter  vector  n  ■  ( n  , 


(yN0,+  ’  °'  SO  ’  ° 

lyH0,+  ‘  ^  SO  ■  ° 

+  .T 

n  )  .  For  any  constant 


K  <  C 


23 


-24 


(where  c23 

(6.45) 


la  defined  aa  in  La  mu  (6.36)  thara  axiata  a  K(k)  so  that 
'  (yjj )”l  <  K(r)  Isf  IP(0,i,ic>,  j  -  1(1)q 

•((yjj)*.*^)'  <  «<*>  «n"tain{e,P(0,i-1,ie}.  j  -  0())q-1 


(6.46) 


,<yij,+l  *  K(“>  ,’1+n,<i',,-1',e)»  J  "  I-1 

,,{yij)  '*ij)'  <  KU>  ,n^,*Ante»p<i^1»M“‘,»'e> J  "  1(,)<> 


i  *h  -1 

p(k,i,x)  -no*  —■) 

r-k  * 


holds  independent  of  e. 

Proofs  Ns  define  u^  -  (yj^)  P(0.1#«>  “  **jP(0,i,]  \  j  -  1(1)q, 


(u^)*  -  (y*j)+P(0,i-1,x)-1,  j  -  0(1)q-1,  i  -  0(t)N-t.  Then  and  vh  satisfy 


'ij 

(6.47) 


“iA  "  Sij<(Ah1,"“h  4  +  Xo 


(6.48)  <DljUJ  -  «:j((<1>X  >  OXV  +  +  ^>"Xq 


(6.49) 


(6.50) 


“iA  “  "iO  4  XjXVh  4  AhXV 


u00  ■  n  '  "NO  “  °'  V80  "  0 

diJ  ■  (1  +  ~>_1'  3  "  H1)q-1,  d’0  -  1,  d^  -  1.  j  -  K1)q 
4*0  -  1  *  0i3  -  diag((dJJ),..#(d’jf\  (d^J)_1,..,(d^j)"1) 


(To  derive  (6.48)  the  equations  u>*  -  «  .  -  w  ,  .  .  and  -  w*  -  at4 

jk  q,q-k  q-J.q-k’  Jk  q.q-k  q-J,q-k 

have  to  be  used.)  Since  ve  can  look  at  (6.49)  as  a  perturbation  of  the  explicit  Euler- 
schesM  the  terminal  value  problee  for  (6.49)  is  stable.  Thus  we  obtain  (using  a 


perturbation  argument  for  (6.47)  and  (6.48))  that  Ks^.v^)!  <  const  Ini  holds  if  k  Is 


less  c2j  (where  e23  is  defned  as  in 


(6.36)).  Because  of 


(6.28) 


-25- 


“  <J<*)  hold*  and  therefore  *(y^)  +  *h  ia  »lao  of  ordar  0(e).  Tharafora  (6.4S) 
holda.  (6.46)  la  proved  analogoualy. 

Proof  of  Theorem  ( 3 . 1 0 ) »  To  complete  the  conatruction  of  the  general  solution  of 
Leh<yh'Xh)  ”  *h  "*  *irst  d«f in®  a  aolution  (yh#«h)  of  the  inhomogeneous  problem  which 
depend*  on  the  »- dimensional  free  parameter  vector  hi 

(6.52)  WVV  "  V  <?00r  *  °'  (yiiO)  +  "  °*  B*(V“h)  "  *  • 

Let  (y^.rj).  (y^,r^) ,  n  ,  n+  be  defined  aa  in  Lemaa  (5.45).  Then  lyh+yh'4Yh'“h+*h**h> 
la  the  general  aolution  of  L^ly^a^)  “  6h*  To  prove  the  exlatence  of  a  aolution  of  0.6) 
and  the  atabllity  of  the  aaaociatad  operator  we  ahow  that  the  equation,  which  ariaea  when 
we  inaert  the  general  aolution  In  the  boundary  condition*  (3.4)  ia  aoluble.  We  obtain 


(f]  \  ^ 

n+'  *  Bi2(yoo,yKO'“oo'*MO)  +  o(ein",n+»  -  B1 


(6.53) 


n  ♦  0(ei(n  ,n  )i)  «  b2  . 


Thia  system  ia  Invertible  and  the  now  of  the  inverae  ia  bounded  uniformly  in  e.  So  there 
exlate  a  constant  cJ4  such  that  l(n”,n+.n)l  <  c24IBI  holds.  Therefore  ,<yh'*h)  1  < 
const (IBI  ♦  ••h,h>  hold*  independent  of  e  and  the  grid  T„. 

(6.54)  Proof  of  Theorem  (4. 1 5 )  1  we  start  with  constructing  an  0(e)  approximation  of  the 
aolution  of  (3.1 1*  Let  (?h«*h>»  yh.  and  (y°,s°)  be  defined  by 

(6.55)  r*(?h.;h)  -  0.  b*(?h,Ih)  -  0,  ®hy‘;  -  8"(f;(tH,2h)^"> 


*I>h  y*  -  8^(f*(th,*h)yJ|) 


(J00\  -1 

]  -  -V  <*00'*ll0>(,,1<*00#S0)<b12<*00#*ll0)  *  b2**00'“n0#0'0'0** 

ywo  ywo 

0  ~  *  0  ~ 
yh  ■  yh  *  V  *h  ■  “h  • 

Mere  V  and  W  are  defined  as  in  (2.7).  Because  of  (2.4)  there  exists  a  aolution  of 

(5.55).  Let  Ujj  and  vh  be  defined  by  ^  “  y^-y*.  v|,  “  *h“*h*  T*,,n  **h  an^  vh 


t 


satisfy 


“h%  “  V»f ( 4  (W 
"  8h,S~<%,\ 4  »r  ‘  w 


4  VW  4  0(€W 
4  82^tth,Yh^  4  0(e,*hV 


“li-h  ‘  8!j‘£‘VUh  4  fr  W  4  VW  4  0(e*ij') 

&  «%»£)  ♦  VW 4  °<e,(W  V  - 0  • 

n  n  n 

*S,  4,#8tM  ‘VW0,  VW  mi(VYh)  *  VVYh)l 

<  c  mx  l(u^,vj|)l  *  **uh“uh'vh”vh'**  (0,0 )  ■  0,  i  ■  1(1)4,  for  mm  constant  c^j 

independent  of  c.  Bacauaa  of  Theorem  (3.10)  we  can  now  uaa  nonlinear  Picard-iteration 

starting  with  «  (0,0)  to  ahow  tha  axlatanca  of  a  eolation  of  (3.1).  We  obtain 

0  0 

l(yh-yh»«h"*h)<  *  0(c)  a Inca  tha  firat  atop  of  tha  iteration  ia  of  order  0(c).  Keller 

(1975)  ahowa,  that  if  a  nonlinear  operator  haa  a  Lipachita  continuoua  linearization  within 

0  0 

a  aphera  {{y0  0}J  than  tha  operator  ia  atabla  in  a  aphere  Kp  <<Th»*h>>  and  P,  and 

0  h  it  1 

tha  atabllity  conatant  depend  only  on  tha  Lipachita  conatant  and  tha  nor*  of  tha  inwaree  of 

0  0 

tha  linear laation  of  tha  operator  at  (y.  ,a.  ) . 

ki  n 

(C.57)  Proof  of  Theorem  (4.1B)»  Let  y0(t),  a0(t),  u0(|>,  vfl(-~)  ***lned  aa  in  (2.1), 
(2.2)  and  let  (y^ii^li  yh  be  defined  aa  above.  Since  (5.55)  la  the  diacrete  counterpart 
of  the  reduced  problea  and  we  uae  a  method  of  order  (q+1 )  we  obtain 

(6.58)  »<v%>  -  VV”  “  0<,,^1, 

by  using  (2.7)  and  standard  arguments.  Bacauaa  of  Lamma  (5.19)  we  Obtain  for  I 
sufficiently  small 

.  i  oh  .  .  8-1  ch 

(6.59)  ly  I  <  oonat  i(1t  —*)“  ,  J  •  l(Dq,  ly..l  <  const  *  (1  ♦  — -“),  j-0,<j-1. 

13  r-0  *  13  t-i  e 

Bat imating  MyJj.sJj)  -  !yo(ti J ,'*0,tij) > '  **  ,(*h'“h,~<*0,th>**0(th>>  4  •*ijl  4 

lp#(t  ).  ♦  iw0(t  )l  gives  tha  raault. 


(6.60)  Proof  of  Lemma  (3 « 22 ) i  Let  (yft^, £),«(th, e) )  denote  the  eolation  of  (1.1) 


restricted  on  the  grid  Th.  Then  we  have 


(6.61) 


eDijy(th'e)  “  43ey,<tio+hi*,e><1* 


If  we  interpolate  y* (t^g+h^s,  )  at  u  v  *  * ' uq  by  a  polynomial  of  degree  q-1  (which  is 
integrated  exactly  by  )  and  use  the  formula  for  the  remainder  term  of  polynomial 
interpolation  we  get 

(6.62)  l(iy  )“l  <  const  chq  ly(q*1)l.„  _  .  . 

XO  1 

Similarly  we  get  I  Uy  )+l  -  0(ehqlytq+1 1 1..  „  ,),  It*  I  -  (Xh?*1  ls<q+2)  I,  .).  It 

follows  from  the  asymptotic  expansion  (2.1),  (2.2)  that  y 1 q+ 1 ' ,  r*q+2)  satisfy 


(6.63) 


t_, 

e*p( — —)  *  exp(X  — -) 

I  (y  ^  (t,e),  z  **  (t,e))l  <  const (1  ♦ - -  —  .  . . ~) 


(6.64)  Proof  of  Theorem  (4,25)»  Let  y(t,e)  be  defined  by 

(6.65)  ey*  -  A<  t,*Q(t , e) >y,  y  (0,e)  -  (JQ(0),  y+(1,t)  -  v*(0) 

where  t„,  uQ.  «0  are  defined  as  in  (2.1).  As  it  is  easy  to  show  y,  w  and  v  satisfy 

•y  ~>lo,(o  1)  "  °<c,»  *y+  ”  wo*(o  I)  "  now  yh  and  *h  be  defined  as  in  the 

proof  of  Theorem  (4.15).  Since  1*. -l„(t.  , 6) I.  «  (MB*1*1)  holds  and  the  initial  and 

nun  n 

terminal  value  problems  in  (5.55)  for  y.  are  stable  we  obtain 

n 


,yh*  (th'€),h  i  con'tl,V‘o(S'E,,h  + 


(6.66) 


♦  max(ehqly<q*1,l..  ,)  t  i  -  0(1)n-1)  -  0(«) 

1 


if  we  choose  the  mesh  according  to  (3.24).  Sumaarising  (6.65)  and  (6.66)  we  obtain  that 
0  0 

**yh'*h*  ”  (W£>'  *o,th'e,,,h  "  holds,  where  (y0, eg)  is  defined  as  in  Theorem 

(3.18).  Therefore  •<yh**h)  ”  (ylt^.e),  »(th,e))l  «  0(ct6)  holds.  So  for  c  and  6 
sufficiently  small  (ylt^.OsU^.c))  lies  within  the  ball  Kp(yh,sh).  80  we  can  apply 
Theorem  (3.15)  and  obtain  that  Ky^.a^)  -  *<^>>1  i  const  4  holds  independent  of 


-28' 


-  It  - 

Proof  of  Bwgjg  (4.27  )i  Ha  choose  hj  “  ttf*  exp(^“)  as  long  as  h^  <  holds  i.e.  as 

ca  1  _  Xh,  _  t.X 

long  as  t  «  r1  tn(-)  holds.  Ry  setting  h.  ■  — p,  t.  ■  — — -  Stop  1  In  (4.25)  btcoof 
i  a  c  l  qc  x  qf 


Chooaa  h.  m  ue  1  as  long  as  t.  <  tn(~)  holds.  (With  a 

X  1  c  <| 


<) 


Given  eQ  and  a 


asauM  that  hs  would  need  M  points,  i.a.  t^  >  tn  —  holds.  If  we  add  1  point  t^+1  > 
+l#(“-  exp(~)).  Therefore  N+1  is  the  number  of  points  we  would  need  to  solve  the  same 

0  eo 

a 

*0  a 

problem  with  e.  ■  e.e  .  Computing  the  sequence  e  •  c  exp<-  — )  shows  that  for 
•  o  nti  n  c 

•3  -2  -la  n 

a  >  10  it  decreases  frost  s  ■  10  to  e  •  10  in  3  steps. 

To  obtain  estimate  (4.29)  we  state 


h^  >  c«q  exp(ia) 

-  i-1  -  la 

t.  >  e4^  J  exp()a)  “  tS*1  * . - 

1  }-0  ea-1 


Therefore 


holds. 


W  <  tn<  tnt-;)  — +  1) 

<X  CO 
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